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GENERAL AIRfOIL THEORY* 
By H. G. Ziissner 

On the assTimption of infinitely small disturljanfies 
the ati.thor develops a generalized integral ecLuation of 
airfoil theory which is applioahle to any motion and com- 
pressiole fluid. Successive specializations yield vari- 
ous simpler integral equations, such as Possio's, 
BirnTJaum ' s , and Prandtl's integral epilations, as well as 
new ones for the v^ing of infinite span with periodic down- 
wash distribution and for the oscillating v/ing with high 
aspect ratio. Lastly, several solutions and methods for 
solving these integral equations are given. 

INTRODUCTION 



There are a nujuher of airfoil theories which hold 
true in two or three dimensions, are stationary or non- 
stationary, and allov/ or disallow for the compressibility 
of air. All these theories have one thing in common: 
They are, strictly speaking, valid only for infinitely 
small disturbances; hence the airfoil must be assumed as 
infinitely thin 'and with infinitely small deflections from 
a regulating surface , the generating line of which is 
parallel to the direction of flight. Then the regulating 
surface itself., can be approximately considered as this 
place of the wing and the area of discontinuity emanating 
from its trailing edge. Up to the present time, a plane 
has been commonly chosen as a regulating surface, but 
this restriction is not necessary. 

following the temporary interest attaching to the 
vortex theory, the introduction of Prandtl's a^cele.rj|.t ipn 
potential mad^ the old pp|;ential theory applicable to air- 
foil theory. The particular advantage of this method over 
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the vortejc method is thiat the cpM^.l^.§3Ai^Jrli:ii3' of air can 
te taken into account. .By this method the wing is re- 
placed hy an arrangement of ac^oustic "radiators." The 
sole essential riestriction of this theory consists in the 
assumption of moderate fields of sound, that is, small 
interferences. Then the classical wave equation 

. □^^^^^^li^^li.-llL.^o (1) 

is .applicahle for the velocity potential and the sound, 
pressure of a quiescent source di s tr ihut ion . Their solu- 
tions have already heen explored in all directions. Its 
application to the moving airfoil is achieved with the 
aid of the well-known Lorentz transformation, the sole 
invariant jjreing the speed of radiation c, which, in the 
case in qv(!estion, is equal to the velocity of sound. 

The setting-up of the integral equations of the air- 
foil theory is a preliminary task, which is definitely 
achieved T5y the subsequent expositions. But this pre- 
liminary work alone accomplishes little without attacking 
the purely mathematical main prohlem, namely, the solution 
of these integral equations without entering into new 
discussions every time regarding the method of derivation 
and its physical significance. 

2. THE VELOCITY POTENTIAL OF THE ELEMENTARY RADIATOR 



The wave equation 0^ = 0 has, as is known, a very- 
simple solution for a spherical wave that spreads .out ra- 
dially at speed c from its source. The solution reads: 

where r denotes the radius, t the time, and f an - 
arbitrary function. Such, a spherical wave is produced 
by an elementary radiator of zero order, which represents 
a simple point source. The velocity potential of radia- 
tors of higher order follows from (2) by partial deriva- 
tion along any coordinate directions. 

The airfoil which is to manifest a pressure differ- 
ence between its two surfaces is best replaced by a su- 
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perppsition of. radliitors of the first order (so-called 
doulDlets) whose axes are normal- to the airfoil.' With n 
denoting the direction of the normals, the potential of 
a douljlet of the superposition is 



The sound .pressure of the, field of- sound is 

P = " p a? 

hence it satisfies eq\aations [Hp = 0 also hecause p~ 
constant for very small di s turlaances . 

3. THE LOREITTZ TSAUSFOEMATIOK 



The arguments so far have dealt with the radiator 
at rest at infinity. |To reach the pressure field of the 
radiator moving at constant speed v < c past the X 
axis, the Lorentz transformation is resorted to. With 
the prime ' ' indicating the transformed system, the 
transformation formulas for a motion along the negativ.e 
X axis of the quiescent reference system read (refer-' " 
ence l) : 



X"i 



=; y = y'; z = z«; t = . d= — (5) 

-/l - V^/CS Jl _ ^2/c3 



On passing to the coordinate x' = X' + vt' in the 
entrained reference system and posing Mach ' s number 
^. = v/ c , ■e(iuation (5) gives: 



X s 



y=y'; z=z«; t = t » Vi - + — (6) 

Vi-p^ ' cVi - ^2 



According to (s), the change to the moving reference 
system .entails the following' suhstitutioh's : 

r^ = -iLlf + yl2 + z!3 (7) 

1 - 
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c " " L c(l - p^) "* c(l ~ 

Inasmuch as the airfoil is to "be situated on a regu- 
lating surface with generating line parallel to flow direc- 
tion, that is, along x direction, the Lorentz transforma- 
tion does not alter the direction of the normals n* 

Later, the knovrledge of the sound pressure of a mov- 
ing douTilet is necessary. According to equations (2) , (s) , 
and (4), the sound pressure for the quiescent douhlet, the 
sequence of the differentiations being changed, is: 

e af 

°ince f is an arbitrary function, can "be 

replaced by another arbitrary function Y; and constant 

factors such as J 1 - can then be omitted or added. 

After completing the Lorentz transformation, equations 
(?) and (8) give for the moving doublet the pressure 

Vft + ■' ^ -t- (1- sg)(y^+ ' 

p(x,y,z,t) =- 7- T 7z=i:r=zzziiiz=rz==zzzzzzzi 

v;here the transformation prime can then be omitted. 

4. fhe Acceleration Potential 

According to Prandtl (reference 2) a new approach to 
the airfoil theory is afforded by proceeding from the 
field of the acceleration vector 

vjhere w is the velocity vector. 

Assuming f rictionless , homogeneous fluid gives the 
Euler equation: 
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grad p 



= ^ + grad 




(10) 



Heace Ja is also a gradient, and an "acceleration 
potential" q) can "be introduced, so that i = grad <», 
Then integration of (lO) gires 



If the flovj is undisturTsed at infinity, then f(t) = 
const. If, in addition, it is restricted to small dis- 
turbances, then p const, and the acceleration 



follov:s from (ll) . 

On the airfoil, pressure p, and hence manifests 
a jump (or discontinuity). At any point outside the 
lifting surface the value p, and hence cp t is continu- 
ous for V < c. 

The acceleration potential cp can "be ap prox imated, 
in the case of compress ihle fluid, from a douhlet super- 
position. This procedure wsuld "be rigorously executed for 
J^ft-c ompr ess ilJl e fluid, "but for i^Hcompressible fluid the 
pressure amplitudes vjould have to he small. This super- 
position can also he assumed on the regulating surface in- 
stead of on the lifting surface, provided the interferences 
are very small. 

Through the superimposed douhlets, the intensity of 
v/hich at any point corresponds to the magnitude of the 
pressure jump on the lifting surface, the acceleration 
potential 9 is defined and follows from (9) and (l2). 
Then the velocity field must he defined from q> . On the 
previous assumption that the flying speed v is constant 
and that the lifting surface moves along the negative x 
axis, the acceleration vector is 




(11) 



= - £■ + const 



(12) 




+ V 



3x 



(13) 



the terms of the second order heing ignored. 



Posing lb = grad q), w = grad follov/ed hy in- 

tegration (13), gives 
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8t * 9x 



<p » + v ~ (14) . 



Given cp(x,y,z,t), the desired acceleration poten- 
tial finally follovrs from (14) at . 



X 



^(x.y.z.t) = i rcp(x',y.z.t - ^ ; )dx' (15) 

— CO 

One particular advantage of this method of represen- 
tation is that the so-called area of discontinuity behind 
the airfoil plays no part because, though it manifests a 
velocity jump, it has no pressure jump and hence no dis- 
continuity in cp. 

An element d a of the surface covered with doublets 
furnishes, according to (9), (12), and (15), the potential 
proportion 



x-x' x'P viSj^+ljlliiizfilll 
^ ^ r ±. V c(l-P^) .^CA^-i^^^ ax* (16) 

— 08 



X 

An- 

d $ 



The analyzed element d cr lies in the zero point of 
our coordinate system, that is, at x, y, z « 0. 

5. THE OEHERAI. INTEGaAL E;^UATI01I OF AIRPOIL THBOEY 

The location of our elementary radiator is now shift- 
ed to any point of the airfoil with the coordinates |, 
y(n) , zCn). Coordinate | is measured parallel to the 
generating line of the regulating surface; hence along 
the X a-xis. Coordinate n is so chosen that, after de- 
velopment of the regulating surface in a plane,, | and 
form a system of Cartesian coordinates, and the sub- 
stitution 

X y z 

X - I y - y(Ti) z - z(tj) (17) 

must be effected in (16). 



The element of the surface of the airfoil is dcr = 
d|dti. Ordinarily Y will still be a function of t. and 
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T| "because of the variation of doublet intensity on the 
airfoil. The differentiation along the direction of the 
normals can lie divided in ^ 

A = ^ ^ + ± i£ = sin a(Ti) ^ + cos a(ti)' ± (18) 
on oy. dn oz dn °^ . °^ 

where angle ait)) "between the direction of the normals 
and the z axis is a function of the curved coordinate r\. 
Posting (17) and (18) in (16) , followed by integration 
over the airfoil, gives the complete term of the velocity 
potential of the airfoil at: 

/ dx'd|dT|(sin a(n)g^coB a^Ti^ r 



$(x,y,z,t)=~y y J dx«d|dT|(sin a(n)g^coB a(Ti) J: 



! . X«-X+e x'p yx'3+(l-p3)[;(y_y(T,))3+(z-j5(Tl))2] 

• yx'2 + (1 - p2) C(y - y(Ti).)« + (z - z(n))*] 

On the other hand, the pressure jump TT on the air- 
foil is proportional to the intensity of the doublets. 
The constant factors in (9) are precisely so chosen that 

n ( |,Ti,t) = pvV(|,Ti,t) (20) 

This conforms to usual practice: lift, positive upward, 
downwash, pos itive ' downward. Literature at times quotes 
the more abstract? downwash, positive upward, in which 
case the prefix of the right-hand side of (30) must be 
reversed.. 

Equations (19) and (20) represent the most general 
integral equation of the airfoil theory .for small dis- 
turbances that can be use'd for computing the pressure 
jum.p • <|> for a given downwash. Equation (19) represents 
a boundary value problem. Its solution rests on the fact 
that, the downwash -on the airfoil itself is given by the 
type of .motion and form, change of the airfoil in first 
approximation. Assujne that n = n(x,.Ti,t) is a small de- 
flection of the airfoil f.r.om the regulating surface in di- 
rection of the normals. Then. the downwash on the airfoil ■ 
the terms of the second order "being neglected - is 



'(19) 
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w{x,y(n) ,z(ti) ft) = 1^ + V 



(21) 



When solving (19), the Kutta flow-off condition must 
"be met. It is met when Y = 0 on the trailing edge. A 
frequently employed formula, which corresponds to this con- 
dition, is: 



^(l.ti.t) - g(n.t) 



a, 



cot 



e 



+ I, aQ sin n 6 



where 



cos 6 = - Y 



and I = 1(ti), half the wing chord. For special pur- 
poses, such as airfoils with circular contour, for in- 
stance, comhinations of spherical functions of the first 
and second kinJ are used, which also satisfy the Kutta 
condition. 

It may "be mentioned that the vortex method also 
arrives at the quantity "Y defined hy (20). .It indicates 
then the density of the hound vortices. This quantity Y 
has, however, a much more general significance than doub- 
let intensity indication since the vortex concept Is re- 
stricted to incompressilsle fluid, while (19) and (20) ap- 
ply equally to compress ihle fluids. 



6. SPECIAL POEMS OF INTBGHAL EqUATIOH (19) 



The general integral equation {19) can lie specialized 
in several directions. 

a) For solving the houndary value problem , the down- 
wash on the airfoil itself is usedi that is , for points 
of origin with the coordinate y = y(Tij_); z = z{r\^) , 
Since solution (2) of the wave equation upon which our 
integral bases applies only otitside the singular point 
r = 0 in view of the linearization, one may not integrate 
through the doublets because then the integral over x' 
becomes divergent.- Hence the downwash on the airfoil 
itself can be obtained only by a limit transition from 
(19), The method nearest at hand consists in carrying 
out the integration with respect to two surfaces at dis- 
tance +e and -e f rom • the airfoil, in forming- the aver- 
age value and proceeding to e — >• 0. This method, applied 
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to ..(19), leads to compircated formulas*. On the other ■■ 
hand, .the same result, i.s achieved if, after the dif.fer- 
entiations, y = 7(111) and z = z(iii),:: are formerly 
written in (19.) and the divergent integral is then so 

■transformed hy partial integration that Cauchy's; princi- 
pal value can he formed. This change makes the formula 

-longer and, in. general, less comprehensive. . So in the. 
following we confine ourselves in most cases to the 
.st.atement -of the simpler divergent integral while tacit- 
ly presuming the further trans for.matipn in. the Cauohy. 
principal value. .If. the integration witJi respect to . x' 
is made numerically, the divergence can be avoided quite 
simply, .by- either integrating from s - I to ■ -00 or 
from +00 . to X -|, depending upon whether x - |< 0 
or > 0. Tor in any event the downwash must disappear at 
infinity. In the divergent method of writing the down- 
wash on the airfoil ultimately assumes the form: 

x'=x-| 

— = w(x,y(Tii) ,z(tii) ,t) = — J dx«d|dTi 



^sina (Ti)^H-cosa(Ti)^ 



\^ V. v(l-ea) c(l - gs) 

yx'8 + (1 - p2) [(y y(n),)s + (z -'z(ti))2] 



y=y(Tii) 

z=z (tii ) 



h) It can he asstimed that the regulating surface is 
a plane, say, the xy plane, ,for instance. Airplane wings 
are usually flat structures. .Then a ss 0 , y('n). = T> , and 
z(n) = 0, With specialiaztion a) z = 0 also. This as- 
sumption itself effects a substantial simplification in 
the eauation form and has been practically always intro- 
duced in the airfoil theory sO; far^ 

c) Some special assumptions regarding" the type of,, 
time rate of change can be made. The linearity of the 
integral equation affords especially simple forms on the 
assumption of harmonic processes with respect to time; 
hence harmonic ^oscillations ;0f the airfoil. The preemi- 
nent importance of the harmonic solutions of the wave 
equation in. phy s ics is an established- fact. .Once harmonic 
solutions of the integral equation are known, solutions 
for any time rate of change of downwash can be found by 
superposition whereby the r-»v9r2^g# Laplace transformation 
plays a prominent part. Proceeding to the limiting case 
of very small oscillation ahd putting ;.Y = ^(I.ti) results 
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in the stationary form of the integral equation. Encum- 
bered with further restrictions, this special case has 
been treated most so far. 

d) Simplified assumptions regarding the shape of the 
airfoil can be m^de.^ Allowing the span of the airfoil to 
increase to -infinity and assuming the downwash to be in- 
dependent of X] gives the same flow process in all the 
planes ti = const, so that the integration with respect 
to 11 can be effected. The form of the integral equation 
so obtained represents a two-dimensional flow process, 
This form is of particular interest because it is solved 
for any downwash- function under the restriction 3=0 

and for stationary flow also with ^ ^ 0. The significance 
of this solution in the so-called vortex filament theory 
will be explained elsewhere. 

On an airfoil of constant chord and infinite span 
both the downwash distribution periodic in r\ 

vii ,r\,t) = Wp( |,t) exp i M- Ti 

and that independent o f r\ can be taken into considera- 
tion.. This defines a "Y distribution periodic in 
The integration with respect to t\ can be effected in 
various cases. ?rom these harmonic solutions superposi- 
tion affords solutions for any finite airfoil of constant 
chord, in motion. The method is therefore applicable 
only to airfoils having a parallelogram (especially a 
rectangle) as contour and pressures that no flow passes 
around the lateral edges of the parallelogram. 

The ultimate aim of the simplifications is the 
change of the surface integral in (19) to a line integral, 
which naturally is more promising for a solution of the 
integral equation. This aim can also be reached by as- 
suming airfoil contours representing coordinate lines in 
simple systems of curvilinear orthogonal coordinates. On 
airfoils of elliptic plan form the confocal elliptic co- 
ordinates permit integration with" the aid of the Lame 
functions, For circular contour there are the spherical 
functions of the first and second type. Examples' can be 
found in the- reports by Kinner (reference 3) , Kriehes 
(r.ef erience. 4) , and Schade (reference 5) . 

e) It is readily apparent from (19) that the assxuhp- 
tion of incompressible fluid affords a substantial sim- 
plification in the form of the integral equation because 
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p = 6 and c = oo . One' eiception' is the stationary flow, 
where, mpreover, the complicated argument of V, which 
starts with t, disappears. Allowance for the compressi- 
hility then consists merely in the reduction of the x 
coordinate in the ratio 6^, the so-called Prandtl 

contracti.oJi., -which represents the exact counterpart of the 
Lorentz contraction in the special relativity theory. The 
downwash therefore should be computed for an airfoil with 
less ch-ord, which, thus- affords hy given pressure jump a 
smaller downwash and, conversely, by given downwash a 
higher pressure jump than for P = 0. 

In the case of stationary flow it forthwith affords 
solutions ^ ^ 0, as soon as solutions for p = 0 are 
available. 

The" oscillating airfoil for 0< & < 1 has been 
treated in only one report (reference 7) , and then, as 
a two-dimensional problem. All other investigations of 
nonstationary airfoil motions are restricted to the assump- 
tion of incompressible fluid. The reason for this is to 
be found in the fact that the general case P ^ 0 v;as not 
approachable by the vortex method solely employed hereto- 
fore until Prandtl introduced the acceleration potential. 



7. TEE AIEPOIL 01 INFINITE SPAH 



The integration with respect to r\ on an airfoil of 
constant chord and infinite span in (19) is dictated by a 
special assumption for the relation of function Y with 
time T. We put 

V(|, T ) exp ip T 

that is, study the harmonic oscillations of the airfoil. 
This assumption implies no loss in, general validity since 
any other-s can he built up from the h&.rmonic solutions by 
superposition. As the ultimate result a closed solution 
eve'n cah be indicated which is applicable to any downwash 
fu^iction without resorting to harmonic analysis. With the 
specializations a), b) , and c) we obtain from (32) by 
derivation according to z the downwash along the z axis 
at : 
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+ 1 +00 X' = X-| 

. r l-x x« ^/x»3+(l-^^.)((y-Tl)^+2^) 1 

. i e.p l^-L-^-^.^- J ^^^^ 

92^ .. yx«2 + (1 - P^)(.(y -. Tl)2 + z2) 

When the new variable 

3 (1 - P ^)(y - n)^ 
_ — + 1 

is introduced into (2S) , "bearing in mind that the part in- 
tegrals with respect to ri TDotvfeen the limits + co to y 
and y to -c» aust Tae of identical magnitude, the inte- 
gration with respect to n can "be effected with the cylin- 
drical function (reference 6) 

00 

exp(-ixu) -f=M=^ = - ^ eJ'^^x) (24) 
and (23) and (24) give the Possio integral equation 

w(x.z.t) = -^4== f ' f dx'd|Y(|,t) exp - 



-co 



I 

X 



1 - P^J 



^ . (^) px,«- . (1 - P^l^ J ( J 

0 c(i - ^^) ' 



Possio (reference 7) does not give this equation explicitly, 
hut merely its kernel in the complex real form in order to 
compute numerical approximate solutions, EqV^) denotes 
Hankel's cylindrical function of the second type, I half 
the airfoil chord. The origin of the. coordinate system is 
at airfoil center. 

After differentiation and Insertion of z = 0 , 
equation (25) gives the downwash at 
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„(x,0,t) = 4^ ' r Y(|,t) exp 1^- (I - x) 



wherein y 



HO.y) = ^ P y /^exp i u H^^^^O.u!) 

is / . ' lU I 



' — 00 

-CO 



^ P y / exp iuHi^^^p u) ^ f-or y>0 



CO 

in 



P y /~ exp (-iu)Hi^®^0 u) for y < 0 (27) 

The absolute value |u| appears in (27) and farther on 
because the argument h&s the significance of a radius. 
With 

lim X H, (x) = ^ 
x-> 0 

it affords for the function R the special values. 



SO.O) = 1 (28) 
R(0,y) = -y / exp i u ^ (29) 



^ du 

— CO 



In stationary flow . v = 0 and (26) and (28) afford 
w(x.O,t) = r yii) (30) 

. , 2 IT . / • ... X - I 



The root factor represents Prandtl's law, according 
to which the c ompres sibility raises the pressure jump in 

the ratio' - 3®. Integral equation (so) was origi- 

nally solved closed by Munk (reference 8) for P = 0. 

In compressible fluid P = 0 and (26) and (29) give. 
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with another integration variable x" =..x - u, the 
Birnhatun integral equation. 

+ 1 00 

Birnbaum (reference (9)) writes his equation in the 
convergent form following from (Sl) after partial inte- 
gration and forming Cauchy's principal value, wherehy the 
sequence of the integrations can be changed. 

Equations (30) and (3l) could also be deduced from 
(23) direct, by putting v = 0 and c = w before inte- 
grating. 

Birnbaum 's equation was solved closed by Kiissner and 
Schwarz (reference 10) with allowance for Kutta's flow- 
off condition, Tor more convenient representation of the 
solution it is expedient to introduce the new variables 

■x=5-lcos6f t=-lcos^ 

and the parameter 

u; » I uj = — - 

V 

Then the solution, on the premise of harmonic time func- 
tion, reads 

TT 

Y(e,t) = i -[^[l+cos^ + T(-iuj) (1-cds d ) 3 cot I 

o 

+ g ?in e ^.^gi^ l-cos(6+'d) |^(^^^)^^ (32) 

cos 9 - cos i l-cos(6-^)-' 

Ho^^^uJ) + i nj^^im) ^ ^ 

•Sim) = ° y-y r-T (S3) 

The complex function T can be computed with the 
aid of tabulated functions and is given in table 1. 
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This resulthas "be en 'applied Ijy the author (reference 
11) with the aid of the superposition principle to any 
downvash function w(dis)f is the path of the airfoil. 

Let w(^,s) =0 f6r s < 0, Then the- solution reads: 



TT 




+ sind In l-cos(e+^) Bw.(d.8) 
l-oos(6--&) 9s 



(1+C08 d)cot jw(^.8) 



2 cos e - CDS' 




(34) 

oo J 

Ui(s) » ■:rT- r exp (ju s dcu (35) 

2Tri / AO) 

-ico S 

The real function Ui(s) can be represented by 
series or a pure real integral over tabulated functions 
(references 11 and 12). Numerical values are given in 
table 2. It should be noted that (34) holds for variable 
flying spoed too, but (32) for constant flying speed only. 

Hext we shall analyze the yawing airfoil of infinite 
span and constant chord. The angle between leading edge 
and axis y, the so-called angle of yaw, is denoted with 
5. The half airfoil chord T is always measured parallel 
to axis X, Putting c = oo- and substituting 

X - y tan 8 = s ; | - tj tan 6 = a 

gives, from (23>, the downwash - 

w(s,0.t) « +1 ^ 

i-T— n r V(cT,t) exp ^(a--x«) /^"'^f^ (36) 

2tt cob3 6 J ■ J ^ V (s-x«)® 

cr=-l x"=cr 

This result differs from (Sl) only in the constant fac- 
tor cos^S, V/ith given dov/nwash the pressure Jxunp of the 
yawing airfoil of infinite span is therefore smaller by the 
factor cosS6 than that of the normally flying airfoil of 
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equal chord measured in direction of flight. This result 
is synonymous with that, in that only the normal component 
of the flying speed, i.e., that located in the plane nor- 
mal to the airfoil, is aerodynamically effective. It is 
readily apparent that in frictionless fluid the tangential 
component, which corresponds to a displacement of the air- 
foil of infinite span along its generating line, must "be 
devoid of aerodynamic effect, and therefore not ascertain- 
ahle by pressure measurements outside the airfoil. 

But, if the pressure is measured in a reference sys- 
tem that is solidly connected with the transversely moving 
airfoil, the Lorentz transformation is dictated if the 
fluid is compres sihle , With Vg denoting the transverse 
velocity measured contrary to direction y, it is neces- 
sary, according to (5) to substitute 

t = + ^gy yp^ 

A - V2*/C^ 

in the integral (26). 

ITezt comes the case of the airfoil with periodic V 
distribution across the span. Let 

'^(e.n.t) = Y(e,t) exp i p. n (37) 

The study is restricted to the special case c = «», 
that is, in compressible fluid. Putting 

r\ - 7 
U s! _■' rrra; 

while noting that the part integrals from +» . to y and 
from y to -o» must be of the same magnitude, affords 
the determinate integral (reference 13) 

CO 

exp(-xu) , = ^ IZU) - Wo(x)] (38) 

Vu^ +12° 

where N is Neumann's cylindrical function 
Struve function. Abbreviating 

Go(x) = Zo(i x) - ^^(i x) 

Gi(x) = Zi(i x) - Ki(i x) - I 




and Z the 



(39) 
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w(x,y,z,t)=i exp iM-y / / dx|cL| V.^i.t + 



— 00 



Since . . 



(40) 



d &o(x) 
dx 



= - 1 



Gi(x) 



the differential follow.ed by insertion of z = 0, x" = 
X - x'' gives the downwash at 

. +1 +00 

w(x,y,0,t) = ^ exp i|j,y J' J y (^|,t 

x'"=| 

E - x" \ , , , , dx" d| , , 

+ i )Gi(-M. X - X" ) : ^— (41) 

V / • )x - x" I 

Integral equation (41) is similar in structure to 
(26) . Its further integration requires a special assump- 
tion for the time function. If harmonic oscillations are 
-involved, a real formula must: be use.d or (4l) must "be. 
written. Since. 

lim' xGi(x) = - (42) 

X-$>-' 0 .TT • • 

the downwash for |j. — ->0 follows from (41) and (42) at 

+ 1 +CO 

' , i=-l x"=|' 
8. THE AIRFOIL WITH HIGH ASPECT RATIO 



According to (23) the downwash of a "flat oscillating 
airfoil in incompressible fluid is 
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w(x,y.z.t) ^ I r r dx'dUnVCt.n.t) exp (x«-x+|) 




i3 



X ^ , (44) 
9z* Vx ' ^+ (y-Ti) ^+2® 

The Introduction of the new integration variable 

X' = - uJiy -n)^ + (45) 

in (44) affords 

w(x.y,z,t) = dldTi-Yd.n.t) exp ^(-x+l) 



''V(y-Ti)2+z2 

On an airfoil with high aspect ratio we find for most 
points of the surface 

X - y - (47) 



Only in the closer proximity of the starting point, 
which itself may be situated on the airfoil, does inequa- 
tion (47) fail to hold unconditionally. 3?iie more slender 
the contour the smaller the error introduced when the lower 
limit of the integral is approximately put at u = 0 in 
(46). following this,the integration can he carried out 
according to (ss) with known functions, ultimately yielding, 
after differentiation, the downwash 



w(x.y.O,t) rr d|dn V(|,T^,t) exp iJi (-x+|) 



(48) 



Equation (48) is substantially simpler than the orig- 
inal (44), Admittedly, it holds only for special downwash 
distributions in x direction; namely, 
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w(x,y,0,t) = w.(y) exp iVj^ - ■^^ (49) 

as is readily seen from a comparison of (48) and (49). It 
affords the downwash function 

This equation can "be solved implicitly only with a 
supplementary assumption for function V(|,n). As the- 
next assumption it is established that the Y distribu- 
tion as a function of i, up to c factor c(ti) depend- 
ent on Ti , shall be such as if two-dimensional flow ex- 
isted in every section. This Is the second essential ap- 
proximate assTimption. The downwash being given by (49) , 
the desired V distribution can be computed from (32) 
and the integral 

iiC-n) 

K(n) = f V(|,Ti) expi-^dl (51) 

13(11) 

formed. Entering (49) in (32) affords with 

i= io^"^^ - COS t3;iju(n) =vI(ti)/v 

V(^.Ti) = C(ti) w (n) exp (- ^ |Q(n)) cot i CJo(ar) 

+ i J^iTu) + (Jo(tlj) - i T ((U)] (52) 

Iq(ti) is the coordinate of the airfoil medial line. For 
airfoil contours that are symmetrical to the n axis, 
to' 0, 

J(j and Ji are Bessel's functions; T is given by 
(33) as quotient of Hankel's functions. Entering (62) in 
(51) , followed by integration with respect to |, gives 

k(t)) = c(ti) ti wCn). Kn) tmr\)) 

f(uJ) = Jo^(<jj)+Ji^(w)+(Jo('ju) - i Ji(w))^ T(uJ) 

If the flow were actually plane, D{r\) would be equal to 1. 

iSText (51) is entered in (50). Since (50), for rea- 
sons outlined previously, is divergent, partial integration 



(53) 



so 
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is necessary for forming Gauchy's principal value. Then 

S(x) . ^ z / o.Cu) ^ (54) 



w(y)-«r(S.K) + ^ . A(^|y - t,|) 



dE(Ti) dT^ 



dt] y - 11 



(55) 



-•b 



By ■unlimited increase of aspect ratio, dK/ dri tends 
toward zero, Ihen exact plane flow prevails and 

w(y) = r(S.K) = \^^]^, <56) 

•nl(y) f(uj(.y)) 

according to (53) and (55), The final form of integral 
equation for the oscillating airfoil of liarge aspect ratio 
follows from (55) and (56) at 

+h 

/^" E (y) 1 Tc/^. ,\ iT^M dT| 

w(y) = r + — / s (— ky - -n ) (5?) 

Til(y)f (w (y) ) 4iT \v'' d n y-Tj 

-h 

Putting u a 0, 

S(0) = 1; f(0) « 2 

gives, from equation (57), the stationary form 

■ +h 

^(y)«.EM_^i rmJll.AJL. (58) 
2TTl(y) 4Tr J d 11 y - Ti 
-h 

This is the well-known Prandtl integral equation of 
"vortex filament theory. " One usually suhstitutes in the 
first term of (58) for tt the measured stationary circula- 
tion constant of a real airfoil section Ci<tt, although 
the whole theory dealt with here applies only to airfoils. 
(Quantity K .in (58) has a very elementary significance; 
namely, 

KM = y Y(i.Ti)di 

that is, the 6o«-cal led circulation of the wing. The sig- 
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nif icance ,of_^^ E . in . ( 57L.J.Ei .less . elementary , Cicala (ref- 
'erence 14). and others have attempted to extend the vortex 
filament theory hy means of vortex concepts to the case 
of the oscillating airfoil and actually arrived at an 
equation agreeahle with (57).. But Cicala arbitrarily 
identifies K with the total circulation (sum of hound 
and free vortices along win,g chord) and accordingly finds 

. = 5 u3 exp i u5 CHo^^^ui) - i E^^^\o5)2 

As concerns the permissihle downwash distribution 
in X direction no data are given at all. So, in this in- 
stance, it is necessary to state, the vortex method leads 
to inaccurate results. Equation (58) cannot he general- 
ized to (57) by vortex theory, although (57) can be spe- 
cialized in (58). 

Equation (57) first gives quantity K(y). But to 
compute the lift requires the bound circulation 

■r(Ti) = f V(e.Ti) d| (59) 

u 

Putting (52) in (59), followed by integration with respect 
to I yields 

r(ti) = C{r\) ttw(ti)1(ti) exp ^ Iq ^^o^^) 

+ i Ji(IIJ) + (Jo(IO) - i Ji(UJ))T(IIj)] (60) 

Comparison of (53) and (60) giyes the bound circula- 
tion 

■E(y) exp (- ^lo(y)) 

^^^^ ^ Jo(5J(y)) - i J^(My))" . ^^^^ 

It is again pointed out that (57) obtained from (44) 
by two approximate assumptions applies only to downwash 
functions of the form (49) , where w(y) can be a complex 
function. Other downwash functions cannot be treated by 
(57) unless. one is satisfied to subtract, a downwash func- 
tion w in the form of (49) from the given downwash 
function Wg and . to integrate the. smallest possible dif- 
ference Wg - w on the assumption of twp-dimensional 
flow, bj means of (32). This method has been used until 
the present time in the practical application of the sta- 
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tionary form (58) without being at all clear Jjhat it in- 
cludes a supplementary assumption going teyond the vor- 
tex filament theory. 

When an attempt is made to apply (57) under the 
cited assumptions approximately to any other downwash 
distrihution w(d,y)» the question of what mean value 
to employ for Tir(y) in. (57) arises. By assuming plane 
flow quantity V can be computed with (32) for any func- 
tions w(^ ,y) . Equation (49) is then so chosen that it 
affords the same characteristic quantity K on the as- 
sumption of plane flow as the given downwash w(^,y) in 
plane flow. Then Vr 

' n' 

K^(y) = l(y)exp^lo(y) / V (8 ,y) exp(-iiDcose ) sine d0 

'Jo 

. V TT TT ^ 

= exp^|o.(y)-~~- yr/'{ci+cos^+(l-cosd)T(u))]cot| 

•Jo 'Jq 

+ LiiB hiHisin^ ini2££llijL4.\w(^,y)xexp(-iiucos e)sin9aad^ 

cosB-cos 8 l-co&(6 - •S) J 

= l(y)exp^-^lo(y)) ( 1-cos d) w( ^,y) d d[Jo ( a))+i J x(ur) 

'^o + C^oi^) - I:Jj;(u3)) Kui)-] (62) 

On the other hand, equation. (49) gives, according to 

(53) 

K3(y) « •tTl(y)w(y) i J^i^Y - i (uT)] 

l^oG) + i Ji(uI) + (Jo("^) - i Ji(UJ)) T(UJ)] (63) 
Putting Ki(y) = K3(y), equations (62) and (63) give 

— (1-cos d)w(^ ,y)dd. 

¥(y) = ° • — exp ^ |o(y) (64) 

Jq(uj) - i Ji(uJ) , 

The numerator of this fraction -represents the known 
integral term Wq that ds decisive for the circulation 
in plane flow,' By translational and rotary oBOillatlbns 
of a f lat ■ airf o il» Wq is equal to the downwash at 3/4 
wing Ohord (rear neutral point)". Por stationary flow the 
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denominat or of -the- f-r act I'on (64) is equal ' t o 1; then w = 
Wq . But this does not hold for nonstationary flow. Put- 
ting, as a check, (49) in (64) identically satisfies (64), 

Allowing for (53), (61), and (64), ecLuation (57) can 
he written in the form 




(l-cosd)w(^,y)d^ = — ; 

0 TTl(y) 

J« - iJ, 



r(y) 



^ Jo -iJi 

+ h 



H, 



Ti7 



+ i H, 



Try 



TIT 



Ti7 



J 



d 

dTi . 



r(Ti) (Jo - iJi) exp 1^ L(ti) 



dri 



y - n 



(65) 



where J and H are functions of tuCy) outside the in- 
tegral and of rD(ri) inside of the integral; r is the 
customary hound circulation. In consequence of the ap- 
proximate assumption regarding the mean value of w, equa- 
tion (65) no longer gives the circulation of the plane 
flov/ exact for infinite aspect ratio, which malces it nees- 
sary to put the first quotient of the first brace equal 
to one. This is apparently as justified as in equation 
(56). 

•Visualizing a stationary periodic downwash field, 
while discounting for the time being the effect of the 
vortices already contained in this velocity field on the 
airfoil, a flat airfoil moving rectilinearly and at con-* 
stant speed through this downwash field, is continually 
subject to a downwash of the form (.49) , whence (57) is 
•fbrthwith applicable. Any stationary d'ovrnwash fields can 
be presented by linear superposition of several such 
fields of different periods. Equation (57) is therefore 
especially suitable for treating gust stresses of airfoils 
with large aspect ratio. 

The complex function S(x) occurring in the kernel 
of the integral (57) is given by (39) and (54). Numeri- 
cal" values- are given in table 3. Because of 



1 im X' 
X — > a> 



*Probably an error... 



(66) 



S4 
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' lim X S(x) =. 1^ ■ 

follows from (54) and (66). 

In rapid oscillations, that is, large values of 
~ |y - t\\ » the induction effect given "by the second term 

of (57) is therefore very small, and in the extreme case 
of high frequency an almost plane flow must be counted on 
even on an airfoil of finite span, if the aspect ratio is 
high. This is an important result. 

Theoretically the solution of (57) can he effected in 
the same manner as the much discussed equation (58), al- 
though the addition of the complex function S makes it 
more protracted. The result, which is, moreover, usually 
encumhered "by the assumption (64) in partial cases, con- 
sists on airfoils with large aspect ratio and at the prac- 
tical values of uj in the order of magnitude of 1.0 only 
in a small correction relative to plane flow. 

In case such refinements of the solution are not 
deemed necessary, it is more appropriate to apply an iter- 
ation method which ties in with the exact integral equa- 
tion (44). A somewhat correct approximate solution 
"^^ifflti) for a given downwash w is afforded from (32) 
on the premise of plane flow, Sutering in (44) gives 

a downwash ^ m. The difference w - Wj^ is then en- 

tered again in (32); is computed and this added to 

V, as correction factor, etc. The convergence of this 
method needs to he checked of course from one case to the 
other, although it should he sufficient in general for 
slender airfoils. 

There is nothing to prohibit the application of this 
method to compressible fluid (^ ^ O) , once the general 
solution of (26) in form of an integral representation of 
the type of solution (32) found for 3=0, is available. 



9. SYSTEMS OF SIVEHAL AIRPOILS 



The general form (19) of the integral equation of the 
airfoil theory comprises the possibility that the surface 
integral d|dn can be extended over severalspatially sepa- 
rated regions of the airfoil. A case in point is the 
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"biplane or the wing with split.:;fl.ai). Of course, Zutta's 
flow-off condition must "blT satisfied for the trailing 
edges of each part of the airfoil. 

Griven a general solution Y = Ij(w) 'of, the integral 
equation for a single airfoil, the problem of two airfoils 
can be splved by succfese ive.' approximation according .fro the 
superposition principle'; .Assume the. given downwash of the 
airfoil (l) aa' Wg^^^ and." of airfoil (2) as . WgA^ = 0, 
?irst' compute- as solution' in fir'st approximation- ' 

As a result of this pressure distribution, there is in 
.conformity with the general integral equation, an aiwajii"s . 
computable downwash field 

which induces oh airfoil. (s) the pressure distribution 
Which in turn creates the downwash field 

..which induce? on.:.'airf oil (l) the. ipressure distribution 

A'y(^) = L^^\aw('>) 
■ -.'.i-J : • ■ ■ . ■ ■ • 

.In the second ■ approximation- the 'pressure distribution of 
airfoil (l) is therefore .given by; V^i' + iiV'^^; and 
that of airfoil (s) by AV^^'. This method of iteration 
can then :be continued. . In -casss where the distance of 
the airfoils is not very small in relation to' the wing 
chord, it converges very well, as proved by Kleinwachter 
(reference 15.) for , stationary flow. 

. The-'iteration method is chiefly opportune when the 
solution Ii_ in continuous form is known. Thus for the 
present it is restricted to two-dimensional flow" with 
c = 9?, for w.hich the .solutions -( 32) and (34) , respec- 
tively,, are available. 



Translation by J. Tanier, 
National Advisory Committee 
for Aeronautics, 
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Table 1. Function r(a»). 



Table H. Function Ui (s). 
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